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SUMMARY 


An experimental investigation has been conducted to determine the effects of rough- 
ness size on the position of boundary -layer transition and on the aerodynamic character- 
istics of a 55° swept-delta-wing model. The tests were conducted at free-stream Mach 

6 7 

numbers from 1.50 to 4.63, Reynolds numbers per meter from 3.3 x 10 to 1.6 x 10 , 
angles of attack from -8° to 16°, and roughness sizes ranging from 0.027-cm sand grit 
to 0.127-cm-high cylinders. 

Near 0° angle of attack, effective roughness Reynolds numbers for the model used 
in this investigation were generally less than existing flat-plate results at similar free- 
stream conditions. Within the frequency range of the thin-film gauge instrumentation 
used in the tests («15 kHz), roughness size had little or no effect on turbulence level at 
the model surface once fully turbulent flow was established. 

Increasing the angle of attack caused the transition to move downstream on the 
model leeward surface and upstream on the windward surface. At a Mach number of 
4.63 and a 16° angle of attack, large regions of laminar flow occurred on the model lee- 
ward surface for all roughness sizes tested. 

An empirical method was derived for predicting the drag of roughness elements 
used in boundary- layer trips and results from the method are in good agreement with 
measured values near zero lift. Experimental results from the tests performed in this 
investigation indicate drag increment due to roughness remains approximately constant 
with changes in angle of attack up to the maximum angles tested. 

For the full range of test variables, grit size had little or no effect on lift coeffi- 
cient and pitching moment. Significant effects of grit size were observed on drag coeffi- 
cient and lift-drag ratio. 


INTRODUCTION 

Limitations in size and unit Reynolds number capability of most continuous -flow, 
supersonic wind tunnels generally require the use of artificial roughness, commonly 
referred to as boundary-layer trips, to promote fully turbulent boimdary-layer conditions 
over surfaces of models. The experimental aerodynamicist is frequently confronted with 
the problem of selecting a boundary -layer trip that is sufficiently large to insure that the 
boundary layer is effectively tripped but not so large that the boundary- layer and local 


flow conditions will be distorted. Selection of grit size is further complicated by the 
current policy of operating wind tunnels at reduced unit Reynolds numbers to meet the 
national goals of reduced energy consumption. 

Numerous investigations have been conducted, with results available in the litera- 
ture, to serve as guides for selecting the correct boimdary -layer trip for given test 
conditions. For example, references 1, 2, and 3 summarize techniques for selecting 
trip sizes for subsonic, supersonic, and hjq)ersonic speed ranges, respectively. Data 
presented in these references as well as most other data in the literature are generally 
limited to zero-pressure-gradient bodies such as cones and flat plates. Very few data 
are available to indicate the effects of pressure gradient, angle of attack, wii^ sweep, 
etc., on roughness required to fix transition. Since most practical wind-tunnel models 
incorporate all these variables, the magnitude of their effects must be determined before 
the zero -pressure -gradient results can be confidently applied. This investigation was 
initiated to provide a data base for some of these effects. An additional objective of the 
investigation was to evaluate the technique of predicting drag of the boundary-layer trips 
described in references 4 and 5 when applied to a three-dimensional body. 

An existing 55° swept -delta -wing model with a 4-percent-thick circular-arc airfoil 
was used in the investigation. This model was one of a family of models previously 
tested in the Langley Unitary Plan wind tunnel (ref. 6). For the present tests which were 
conducted in the same facility, boundary -layer transition location and force and moment 
data were obtained for roughness sizes ranging from 0.027-cm sand grit to 0.127-cm- 
high cylinders at Mach numbers from 1.50 to 4.63, Reynolds numbers per meter from 
3.3 X 10® to 1.6 X lo"^, and angles of attack from approximately -8° to 16°. Both surface 
thin-film gauges (ref. 7) and the sublimation technique (ref. 8) were used to define the 
location of transition. The force and moment data were obtained by using a standard 
six -component strain-gauge balance. 


SYMBOLS 

projected frontal area of individual roughness element, m^ 
Ca model axial -force coefficient, 

Cp, model drag-force coefficient, Q— ^ 

>oo ^ 

Cd c drag coefficient of element used in transition strip, 
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g(A,Mj 
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L/D 
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N 

R 

R’k 

^k,eff 


drag coefficient of isolated cylinder or isolated roughness element, 

Drag force 

^coAc 

model drag coefficient at zero lift, 
model lift coefficient, 

pitching-moment coefficient. Pitching mom _ent 

^ooSc 

mean aerodynamic chord, 36.35 cm 

root mean square of voltage fluctuation, mV 

term in equations (4) and (5) to account for sweep angle and Mach number 
effects 

roughness height, cm 
lift -drag ratio 

length of transition strip, cm 

free-stream Mach number 

number of roughness elements 

free-stream dynamic pressure, N/m2 

free-stream Reynolds number per meter, \5/v, m“t 

Uk 

roughness Resmolds number, 0.01 — 
roughness parameter, equation (1) 

roughness Reynolds number required to move transition near vicinity of 
boundary -layer trip, 0.01 ^ 
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Reynolds number based on chordwise distance from leading edge to location 
of boundary-layer trip, 0.01 

Reynolds number based on chordwise distance from leading edge to location 

Uxf Ux+ 0 

of transition, 0.01 —rr~, 0.01 — 

^ V ^ i; 

2 

model planform area, 0.2045 m 

roughness spacing measured from center to center of elements, cm 

adiabatic wall temperature, K 

temperature at roughness height, K 

wall temperature, K 

free -stream velocity, m/sec 

width or diameter of roughness elements, cm 

chordwise distance from leading edge, cm 

location of transition strip, cm 

location of roughness induced transition (peak root-mean -square voltage 
from thin-film gauge), cm 

location of natural transition (peak root -mean -square voltage from thin -f ilm 
gauge), cm 

angle of attack, deg 

boundary -layer thickness at x = xjj, cm 

boundary -layer displacement thickness at x = xj^, cm 

value of R'jj where x^ = xjj 
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A 


leading-edge sweep angle, deg 


u free-stream kinematic viscosity, m^/sec 

w exponent in viscosity-temperature relation, equation (1) 


APPARATUS AND METHODS 
Tunnel 

The tests were conducted in both the low and high Mach number test sections of 
the Langley Unitary Plan wind tunnel which is a variable pressure continuous flow facil- 
ity (ref. 9). The asymmetric sliding-block nozzles lead to the test sections and permit a 
continuous variation in Mach number from about 1.5 to 2.9 in the low Mach number test 
section and from about 2.3 to 4.7 in the high Mach number test section. 

Model 

Details of the model used in the investigation are shown in figure 1. The model 
consists of a clipped 55° swept delta wing with a 4 -percent -thick circular-arc airfoil. 

A body of revolution was added symmetrically about the wing center line to provide a 
housing for the strain-gauge balance. The base diameter was 5.08 cm and was the 
minimum diameter required to house the balance. 

TEST CXINDITIONS 


Tests were conducted at the following test conditions: 


Mach 

number 

Stagnation 

temperature, 

K 

Stagnation 
pressure range, 
kN/m2 

Reynolds number 
range, 
m"l 

1,50 

339 

27 to 133 

3.3 X 10® to 1,6 X 10^ 

2.36 

339 

38 to 189 



2.86 

339 

49 to 246 



3.95 

352 

92 to 460 



4.63 

352 

126 to 631 


r 


The dewpoint temperature at stagnation conditions was maintained below 239 K to insure 
negligible condensation effects. 
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The effectiveness of transition strips containing roughness sizes ranging from 
0.027-cm sand grit elements to 0.127-cm high cylinder was investigated. The size of 
the sand elements was determined by sifting sand elements of various sizes through a 
stack of U.S. Standard sieves. The sieves were arranged according to size, with the 
largest-mesh sieve at the top of the stack and the sieves with progressively smaller 
mesh toward the bottom. The size of the sand elements resting on a collecting sieve 
was determined by averaging the size of the mesh openings of that collecting sieve and 
the size of the mesh openings of the preceding sieve through which the elements had 
passed. The sieve number, mesh openings, and nominal particle sizes are given in the 
following table: 


Sieve 

number 


Mesh 

opening, 

cm 

Nominal particle 
size collected, 
k, cm 

20 


0.0841 



25 


.0711 

0.0776 ± 

0.0065 

30 


.0589 

.0650 ± 

.0061 

35 


.0500 

.0545 ± 

.0045 

40 


.0419 

.0460 ± 

.0041 

45 


.0351 

.0385 ± 

.0034 

50 


.0297 

.0324 ± 

.0027 

60 


.0249 

.0273 ± 

.0024 


The cylindrical elements were machined from small copper rods and were 
0.127 cm in diameter and 0.127 cm in length. All transition strips were located 1.02 cm 
aft of the leading edge measured in a streamwise direction. For the transition strips 
containing the smaller roughness elements (k = 0.027, 0.032 cm), the elements were 
randomly spaced in a band (approximately 0.159 cm wide) parallel to the leading edge. 
All other transition strips contained single spaced elements that were arranged in a line 
parallel to the leading edge and located at intervals determined from the equation 


cos A 
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MEASUREMENTS AND INSTRUMENTATION 


Strain-Gauge Balance 

Aerodynamic forces and moments were measured by means of a six-component 
electrical strain-gauge balance housed within the model. The balance was rigidly 
fastened to a sting support system. Balance -chamber pressure was determined by 
averaging measurements from two static pressure orifices located in the vicinity of the 
balance. 


Thin -Film Gauges 

After completion of the force and moment tests, the model was instrumented with 
surface thin-film gauges. Each gauge consisted of a thin platinum coating on the end of 
a quartz rod. At 273 K the resistance of the gauges was approximately 5.5 ohms. 

Gauges were mounted at 10 locations on the model by cementing the quartz rods in 
0.153-cm-diameter holes that were drilled perpendicular to the model surface (fig. 1). 
That end of the quartz rod with the platinum film was located flush with the model upper 
surface. 

The gauges were operated in a constant temperature mode at an overheat ratio of 
approximately 1.5. A switch assembly was used to connect each gauge into a single- 
channel anemometer system, and the bridge voltage output fluctuation of the system was 
measured by using a true root-mean-square (rms) meter. 

Typical rms bridge voltages from the present tests for the thin-film gauges are 
shown in figure 2 for both natural boundary -layer transition and roughness induced 
transition. Oscilloscope traces of the bridge voltage fluctuations are also shown for 
comparison with the measured rms voltages at selected Reynolds numbers. The varia- 
tion in rms voltages shows an increase in fluctuation level with increasing Reynolds 
number from a level corresponding to laminar flow at the lower Reynolds numbers to a 
level corresponding to turbulent flow at the higher Reynolds numbers and a peak value 
corresponding to transitional flow at the intermediate Reynolds numbers. For the 
present tests, transition Reynolds numbers are based on the unit Reynolds numbers 
corresponding to the peak fluctuation level, which as discussed in reference 7, should 
give results consistent with values based on peak surface pitot pressure or peak surface 
temperatures. 

For the case of natural transition, figure 2(a) shows large positive voltage pulses 
on the oscillograph trace at Reynolds numbers slightly less than the peak voltage 
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Reynolds number; large negative pulses are shown on the trace at Reynolds numbers 
slightly greater than the peak voltage Reynolds number. The positive pulses are believed 
to result from turbulent bursts intermittently increasing the heat transfer rate from the 
gauge. These bursts require an increase in bridge voltage to maintain the constant 
temperature condition. Conversely, the negative pulses are believed to result from 
intermittent laminar flow occurring at the gauge in an otherwise predominantly turbulent 
boundary layer. With further increases in Re 3 molds number, the location of transition 
moves upstream of the gauge. Random positive and negative fluctuations representative 
of turbulent flow then result. 

For the case of forced (roughness -induced) transition, the peak in rms voltage 
indicative of transition was much lower than the peak for natural transition (fig. 2(b)), 

In some cases no peak at all was observed; however, a definite change in rms voltage 
level indicative of the change from laminar to turbulent flow always occurred when 
transition moved upstream of the gauge. For these cases, the Resmolds number corre- 
sponding to the initial onset of the increased voltage level was defined as the transition 
Reynolds number. Although the measurements of reference 7 always indicated a distinc- 
tive peak for forced transition, the peak occurred just ahead of the initial onset of the 
increased fluctuation associated with fully turbulent flow over a very narrow band of 
Reynolds numbers. Therefore, the definition of transition Reynolds number for the 
present data where no peak was observed is not believed to be significantly inconsistent 
with the technique developed in reference 7, 

Sublimation 

Limited tests were also conducted using the sublimation technique for detecting 
the location of boundary -layer transition. This technique is described in detail in 
reference 8. 


CORRECTIONS 

For the force and moment tests, angles of attack have been corrected for tunnel- 
flow misalignment and for deflection of the sting and balance caused by aerodynamic 
loads. Drag coefficients were adjusted to correspond to free -stream static pressure 
conditions at the base of the model. 

During the thin-film tests, the model was placed at the same location in the test 
sections as it was for the force balance tests. Therefore, the tunnel -flow misalignment 
determined from the force balance tests was used to correct angle of attack for the thin- 
film tests. No corrections were made to angle of attack for aerodynamic loads during 
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the thin -film tests; however, this effect was minimized by the use of a sting that was 
shorter and larger in diameter than the sting used for the force tests. This sting was 
rigidly attached to the model. 


RESULTS AND DISCUSSION 
Thin -Film Data 

As shown in figure 1, the model was instrumented with thin-film gauges at two 
spanwise stations; however, for the sake of clarity and consistency, results are presented 
only for the more extensively instrumented inboard station. When transition was located 
within the instrumentation regions of the outboard station, these results generally agreed 
well with measurements obtained at the inboard station. This agreement indicates that 
the transition front paralleled the wing leading edge. ’ This parallelism was, in fact, 
substantiated by the sublimation data. 

Thin-film data are not presented for Mach 1.5. During the initial thin-film runs 
at this Mach number, it became apparent that the subsonic leading edge permitted 
disturbances created by the forward gauges protruding through the lower surface to 
contaminate the flow over the upper surface. Therefore, further testing of the thin-film 
gauges at this Mach number was discontinued. At the higher test Mach numbers, the 
leading edge was supersonic, and no indication of contaminated flow over the upper 
surface was observed. 

Shown in figure 3 are natural transition results for the range of unit Reynolds 
number and Mach numbers from 2.36 to 4.63. Limited sublimation data that were 
obtained are shown as the primed symbols and generally agree well with the thin-film 
data. The variation of g with unit Reynolds number and Mach number is somewhat 
similar to trends reported in the literature for flat plates and cones. The variation 
of Rx ^ with unit Reynolds munber at the higher Mach numbers, however, is not 
consistent with existing flat -plate and cone data. Extensive wind-tunnel testing of these 
zero pressure gradient bodies has shown that R^ t increases with increasing unit 
Reynolds number contrary to the results shown in figure 3. This unit Reynolds number 
effect is generally attributed to wind-tvmnel noise radiated from the nozzle sidewall 
boundary layer. Several investigators (see, for example, ref. 10) have derived empirical 
equations that correlate transition data from various wind timnels that collectively have 
a wide range of radiated noise levels. Unpublished transition data from a hollow cylin- 
der model and a 10° cone model tested in the Langley Unitary Plan wind tunnel when 
expressed in terms of the correlation parameters in reference 10 are in good agree- 
ment with data from other facilities. Therefore, the trends shown in figure 3 probably 
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do not result from tunnel effects but from three-dimensional effects created either by 
the delta wing itself or by an interaction between the flow field over the delta wing and 
the body section that houses the balance and sting. 

Shown in figure 4 is the effect of angle of attack on the location of natural transition 
for Mach numbers from 2.36 to 4.63. The ordinate parameter is the distance to transi- 
tion at angle of attack nondimensionalized by the distance to transition at 0° angle of 
attack. The results show that on the windward side of the model, transition moves 
upstream toward the wing leading edge with increasing angle of attack; whereas on the 
leeward side, transition moves downstream with increasing angle of attack. The results 
also indicate that the nondimens ional transition location parameter becomes less sensi- 
tive to angle of attack with increasing Mach number. 

Presented in figure 5 is the effect of grit size on the location of bovuidary -layer 
transition. As expected, the effectiveness of a given grit size for inducing transitional 
flow decreases with increasing Mach number. At Mach numbers 2.36 and 2.86, data 
were obtained for k ranging from 0.027 cm to 0.055 cm; at Mach numbers 3.95 
and 4.63, results were obtained for k ranging from 0.055 cm to 0.127 cm. At Mach 
number 2.36 considerable data scatter occurred and a clear trend of the variation of 
with k was not established. At the higher test Mach numbers, the effect of k is 
readily apparent from the data and consists of the anticipated upstream movement of 
transition with increasing k. Sublimation data, shown as the primed symbols in fig- 
ure 5, generally agree well with the thin-film data. 

In reference 11, roughness -induced transition data for flat plates and cones are 
correlated by using a roughness parameter defined as 

0.5+0} 

I (1) 

Since must be calculated at a station Xj^ rather than at x^, the investigators 

reasoned that the effectiveness of the roughness would vary according to the distance 
between roughness and transition. The quantity 


R'k = Rkf^ 


"w 


pT pT/^R'k\ 
^,0 Pt,o^ ^ / 


( 2 ) 


was selected to represent these relative positions for flow having zero pressure gradient. 
The constant e represents the value of R'j^ where x^. = xj^. In references 3 and 12, 
data are correlated using the parameters defined by equations (1) and (2) in the form 
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(3) 



^t,0 


Xk /R'k’' 


where a value of e was selected that best correlated the data of references 3 and 12 
for a given Mach number. Data from the present tests were also correlated using the 
parameters from equation (3). Values of e of 1500 at Mach 2.36 and 2.86 and values 
of e of 2000 at Mach 3.95 and 4.63 seemed to give the least data scatter. These 
results are shown in figure 6 for the test range of Mach numbers for values of k 
greater than the boundary -layer thickness calculated by assuming zero pressure gradient 
conditions over the distance This latter restriction (k ^ 6^^) was imposed due to the 

sensitivity of R'k to k for k < 6k and the unknown error in 6k resulting from the 
assumption of a zero pressure gradient. The reduced values of e used in the correla- 
tion of the present data in addition to the fact that the data level falls below the flat -plate 
data indicate that roughness more effectively trips the boundary layer for the present 
configuration. 


Effective roughness Reynolds numbers for the present tests are shown in figure 7 
and are compared with existing flat -plate data. The effective roughness Reynolds 
numbers for the present tests are defined as the product of roughness height and the 
free-stream unit Reynolds number required to move transition to the first thin-film 
gauge located downstream of the transition strip. Based on the data shown in figure 5, 
this definition should be consistent with effective roughness Reynolds numbers determined 
from the "knee -of -the -curve" technique discussed in reference 3. Results from the 
present tests fall below both flat -plate curves shown in figure 7; however, the curves 
representing data of reference 3 are closer to the present data, particularly at the higher 
Mach numbers. The differences in the values of Rj^ for references 3 and 12 could 
be due to different techniques of defining R^ or possibly to an effect of R^ on 
Rk gjj that is shown in figure 6 of reference 3, but is not apparent from the data of 
reference 12. 


The effect of angle of attack on the location of roughness -induced transition is 
shown in figure 8 for Mach 2.36 to 4.63 at a constant unit Reynolds number per meter 
of 6.56 X 10®. Results obtained on the leeward side of the model are presented in fig- 
ure 8(a) for the minimum and maximum roughness size tested at this Reynolds number 
and at angle of attack. With increasing angle of attack, the results show a downstream 
movement of transition which increases with increasing Mach number and decreasing 
roughness size. The results show that at Mach 4.63 even the 0.127-cm-high cylinders 
were not sufficiently large to retain the location of transition near the vicinity of the 
transition strip for any significant range of angles of attack. This sensitivity of transi- 
tion location with angle of attack could significantly affect the onset of boundary -layer 
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separation that occurs in regions of adverse pressure gradients as, for example, those 
created by aerod)mamic control surfaces. Therefore, caution should be applied in 
interpreting aerodjmamic data (such as control effectiveness, etc.) obtained in wind 
tunnels at high Mach numbers and at large angles of attack. The increase in the extent 
of laminar flow on the leeward surface also results in a reduction in skin-friction drag; 
however, as discussed subsequently, the skin -friction drag at an angle of attack is 
apparently dominated by the friction drag on the windward surface. The effect of 
increasing the angle of attack on the position of transition for the windward side of the 
model (fig. 8(b)) consists of a forward movement toward the transition strip. Therefore, 
if the boundary layer is adequately tripped near the transition strip at a = 0°, only 
small changes in the location of transition should occur on windward surfaces with 
increasing angle of attack. 

Shown in figure 9 is the effect of k on the rms voltages at instrumentation 
locations ranging from the most forward gauge to the most aft gauge. No significance 
is placed on the magnitude of these measurements other than to indicate that once fully 
turbulent flow is established, roughness size has little or no effect on the turbulence 
intensity at the wall. This lack of sensitivity to roughness is somewhat surprising since 
numerous investigators (see, for example, ref. 13) have shown that large lateral shear 
gradients occur downstream of large roughness elements. Part of this lack of sensitivity 
could result from the limited frequency response of the anemometer system used for the 
present tests. The upper frequency limit (-3 dB) of the system as determined by the 
conventional square wave technique was approximately 15 kHz. 

Force Balance Data 

The effects of roughness on drag coefficients at zero lift are shown in figure 10 for 
the test range of Mach number. The shape of the ^ curves are strongly influenced 
by the position of boundary -layer transition. For the smaller grit sizes and lower 
Reynolds numbers at Moo = 2.36, transition is located at some distance downstream from 
the trip location, and reduced values of Cj^ ^ are obtained. As Reynolds number is 
increased, transition moves upstream with turbulent flow covering a larger surface area. 
An increase in friction drag then results and hence an increase in q. As transition 
approaches the region of the trip, a peak occurs in the Cj^ ^ curves and further 
increases in Reynolds number result in a decrease in Cq q. Increasing grit size 
results in a reduction in the peak Cjj ^ Reynolds number and also an increase in the 
level of Cj^ ^ as a result of the grit wave drag. The minimum test Reynolds number 
at Mach 1.5 was not sufficiently low to define clearly a peak in the Cj^ ^ curves for the 
range of grit sizes tested. 
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The crosshatched band crossing the Cj^ ^ curves of figure 10 at g 2,36 
represents Reynolds numbers obtained from figure 5 at which transition is located at 
the first thin -film gauge (x = 2,54 cm) downstream of the transition strip. At all three 
Mach numbers, this crosshatched band crosses the Cj^ curves at Re5molds numbers 
greater than the peak Cj^ ^ Reynolds number, Reynolds numbers required to move 
transition to the second instrumentation station (x = 5,08 cm), as determined from 
figure 5, more closely approximate the peak Cj^ ^ Reynolds numbers of figure 10, 
These results imply that zero-lift drag measurements representative of fully turbulent 
flow can be obtained with transition located at some small distance downstream of the 
transition strip. Conversely, the results indicate that the peak Cj^ ^ Reynolds number 
is not a very good indication of when boundary -layer transition occurs at the transition 
strip. 

It is apparent from the results shown in figures 5 and 10 that at the higher Mach 
numbers large roughness elements are required to induce fully turbulent flow. Under 
these conditions, the extraneous drag created by the elements becomes a significant 
part of the total model drag near zero lift. In the literature (see ref, 2, for example) 
several methods are described to determine the magnitude of this element drag at the 
low supersonic Mach numbers; however, these methods generally require additional 
wind-tunnel testing either through a range of Reynolds numbers or through a range of 
element sizes. Because of the difficulty with fixing transition at higher supersonic 
Mach numbers, and because of the Reynolds number limitation of most supersonic wind 
tunnels, the application of these methods at the higher Mach numbers could result in 
inaccurate results as well as an excessive amount of wind-tunnel operation. 

An empirical method for predicting the drag of roughness elements used in 
boundary-layer trips is described in the appendix. This method is based on drag 
measurements of small cylinders presented in references 4 and 5; its application does 
not require additional wind-tunnel tests. Equations for predicting grit drag increment 
as derived in the appendix are 

ACn = (5xl0-8)|^g(A,M„ 

and 

ACj3^(2.5x 10-^)| kg(A,Mj 
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where 


g(A,M^) = (cos A) Ql - (0.0478 + 0.0430M^)cos A - (0.0341 + 0.0307Moo) cos^ a[] 

[4 '' ^ 


X / cos I /jrt 


, — + 0.0439^1 

Vcos A + 1.992 /J J 


Equations (4) and (5) are derived for cases where k « w, and spacing between grit 
elements is defined by the equation 


s = 


4k 

cos A 


For the present model (A = 55°, S = 0.2045 m^, I = 243.7 cm), equations (4) and (5) 
reduce to the following simplified equations: 


AC 


D 


= (0.00265 - O.OOOlOR^x,) 



and 


ACq = (0.0133 - 0.0005Mjk 



( 6 ) 


(7) 


Calculated grit drag increments from equations (6) and (7) are compared with measure- 
ments from the present tests in figure 11. Since equations (6) and (7) only predict the 
increment in drag coefficient due to the transition strips, it is necessary to establish the 
grit -free turbulent drag level of the model and add to this value the predicted grit drag 
increments before comparisons can be made with the measured total drag values. In 
figure 11, the grit -free turbulent drag levels were approximated by fairing a straight 
line though the data points for k/s ^ 5 to the ordinate axis (k2 = 0 cm^). In general, 
the predicted grit drag increments agree well with measurements for Mach 1.5 to 3.95 
at R = 1.15 X lo'^^ m"^ (fig. 11(a)). Also, the few data points obtained for > 5 

indicate a change in slope in the Cp ^ curves in the vicinity of k/d = 5 as predicted. 

The predicted grit drag increments at = 4.63 and R = 1.15 x lo"^ m“^ are 
considerably less than those measured. For this Mach number and Reynolds number, 

Cr q measurements for the two largest roughness sizes only were used to determine 
the grit -free drag level. Even for these very large grits there is some doubt that the 
k = 0.078 cm grit induced fully turbulent flow at this Reynolds number. A small error 
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in Cj^ Q for the k = 0.078 cm grit results in a large error when extrapolated to 

= 0 cm^ to determine the grit-free drag level. If this lack of fully turbulent flow is 
the reason for poor agreement, then better agreement would be expected at the higher 
Reynolds numbers. To determine whether this is the case, predicted and measured 
grit drag increments are compared in figure 11(b) for R = 1.64 x 10 m~^ and much 
better agreement is obtained. Therefore, the results presented in figure 11 suggest that 
the prediction method gives realistic estimates of grit drag increment for zero -lift 
conditions, at least for the range of test variables of this investigation. 

The data of reference 14 indicate that if fully turbulent flow occurs over a model 
near zero lift, the drag increment due to roughness is relatively insensitive to changes 
in angle of attack or lift. To determine whether this phenomenon occurred for the test 
conditions of the present tests, the data were examined as shown in figure 12. Results 
are presented in the form of in lieu of Cj^ since the drag due to lift of the present 

model results in very large drag coefficients at the larger angles of attack with corre- 
sponding insufficient data resolution to determine the magnitude of the grit drag. Also, 
for the maximum angles of attack of the present tests, and for the roughness 

elements differ by only approximately 5 percent. The results shown in figure 12 indicate 
that the grit drag increment remains approximately constant for the test range of angle 
of attack. Limited thin-film data obtained at R = 1.312 x lo"^ m”^ and = 4.63 
indicate that at the larger angles of attack, large regions of laminar flow occurred on 
the model leeward surface for the smaller roughness elements similar to the results 
shown in figure 8(a). This increase in extent of laminar flow would be expected to reduce 
the skin-friction drag which then results in an increase at the larger angles of attack in 
the difference between values for the smaller roughness elements and the larger 

roughness elements. The constant increment in for the different roughness sizes 

through the range of angles of attack shown in figure 12 indicates this anticipated 
increase did not occur. This lack in sensitivity of the axial-force measurements to the 
increase in extent of laminar flow over the leeward surface possibly results from the fact 
that the skin -friction drag over the model windward surface (flow over which is fully 
turbulent as indicated by the thin-film data) is much larger than that which occurs on the 
leeward surface. The windward turbulent friction drag apparently dominates the total 
friction drag. Preliminary estimates based on two-dimensional considerations indicate 
that the friction drag on the windward surface of the present model at a « 16° and 
Mach 4.63 is over an order of magnitude greater than that which occurs on the leeward 
surface for fully turbulent conditions. 

The effects of grit size on measured aerodynamic characteristics of the present 
model are shown in figure 13 for the minimum and maxirhum Mach numbers of the 
investigation at R = 6.56 x 10® m“^. In general, these results show that grit size had 
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little or no effect on Cl and For angles of attack ranging from 0° to values 

corresponding to maximum L/D, increasing grit size results in an increase in 
and a decrease in L/D, as would be expected. These grit effects result from the 
increase in grit wave drag with increasing k and in some cases to the increase in 
extent of turbulent flow with increasing k. The results shown in figure 13 indicate that 
the consequence of using large roughness elements depends on the particular aerody- 
namic characteristic parameter of interest. For example, if Cj^ and are of 

primary concern, then roughness size is of little importance; whereas if and L/D 

are the principal measurements, then the effects of roughness size must be taken into 
consideration. Also, as previously discussed, the difficulty associated with obtaining a 
fully turbulent condition on the model leeward surface at angles of attack could signifi- 
cantly alter separation characteristics ahead of control surfaces. This effect must be 
considered if measurements such as control effectiveness are of primary interest. The 
model used in the present investigation did not have movable control surfaces and, 
therefore, control effectiveness measurements were not obtained. 

SUMMARY OF RESULTS 

An experimental investigation has been conducted to determine the effects of rough- 
ness size on the position of boundary-layer transition and on the aerodynamic characteris- 
tics of a 55° swept- delta-wing model. The tests were conducted at free-stream Mach 

fi 7 

numbers from 1.5 to 4.63, Reynolds numbers per meter from 3.3 x 10° to 1.6 x 10 , 
angles of attack from -8° to 16°, and roughness sizes ranging from 0.027-cm sand grit 
to 0.127-cm-high cylinders. The results are summarized as follows: 

1. Near 0° angle of attack, effective roughness Reynolds numbers for the model 
used in this investigation were generally less than existing flat -plate data. 

2. Within the frequency range of the thin -film gauge instrumentation used in the 
present test («15 kHz), roughness size had little effect on turbulence level at the model 
surface once fully turbulent flow was established. 

3. Increasing angle of attack resulted in transition moving downstream on the model 
leeward surface and moving upstream on the model windward surface. At Mach 4.63 and 
angle of attack of approximately 16°, large regions of laminar flow occurred on the model 
leeward surface for all roughness sizes tested. 

4. Results from empirical equations, derived from the data of NASA Technical 
Notes D-7369 and D-7812, for predicting the drag of roughness elements used in boundary- 
layer trips are in good agreement with measurements from the present tests. 
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5. Experimental data from the present tests indicate grit drag increment remains 
approximately constant with changes in angle of attack up to the maximum angles tested. 

6. For the full range of test variables, grit size had little or no effect on lift coeffi- 
cient and pitching moment. Significant effects of grit size were observed on drag coeffi- 
cient and lift-drag ratio. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
October 26, 1977 
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APPENDIX 


DRAG PREDICTION OF BOUNDARY-LAYER TRIPS 


Described in this appendix is a method to predict the drag of roughness elements 
used for boundary -layer trips based on the data of references 4 and 5. Although the data 
of references 4 and 5 are for small circular cylinders, it is assumed in this paper that 
the data are equally applicable to sand grit elements. 


In order to simplify the present method, the cylinder results from references 4 
and 5 were converted to the parameter Cp C and were used to empirically 
derive the following equation: 



(Al) 


Data from reference 5 are compared with calculations from equation (Al) in figure 14 
for Mach 2.30 to Mach 4.60 and sweep angles from 0° to 60°. Good agreement is shown 
for this range of variables. 

In order to approximate Cp jc, the results of reference 4 for the case of a 
cylinder having equal height and diameter were approximated by the equations 


^D,IC 



and 


^D,IC 


1.0 



(A2) 



(A3) 


and as shown in figure 15 these approximations are in general agreement with the experi- 
mental data. For the cylindrical elements of the present tests k = w, and although the 
sand elements were of random shapes, it is believed that this equality is also approxi- 
mated for the average element. 
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The roughness element drag coefficient is obtained from equations (Al), (A2), 
and (A3) combined with the equation 


-D,C 




(A4) 


The total drag force for all elements in the transition strip can be approximated from 
the equation 


Drag force = Cd^c^coAcN 


(A5) 


where A^ is the projected frontal area of a typical element and N is the total number 
of elements. 

The increment in model drag coefficient due to roughness drag is calculated from 
the equation 


_ R oughness drag force 
^ ^ qooS 


(A6) 


which when combined with equations (A2), (A3), (A4), and (A5) gives 


AC^ = 0.2N^ 


'D,C 


D 


® 5* Cp) ic 


< 5 


(A7) 


and 


\ ^D,C 




D,IC 


> 5 


(A8) 


where from equation (Al) 


■'D,C 

'D,IC 




Since it was previously assumed that 
k = w 
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then 


Ac = O.OOOlk^ 


(A9) 


Also 

N = I (AlO) 

where I is the total length of the transition strip and s is the spacing between 
elements from center to center. Customarily in many supersonic wind tunnels, grit 
spacing is determined from the equation 

cos A 

which when substituted into equation (AlO) gives 


N = ^ cos A (All) 

Substituting equations (A9) and (All) into equations (A7) and (A8) results in the following 
equations for grit drag increment: 


ACj. = (0.000005) I ^ g(A,Mj 
® ® k 


and 


ACj 3 = (0.000025) I k g(A,M<„) 



(A12) 


(A13) 


where 


g(A,M ) = (cos A) Cl - (0.0478 + 0.0430MJ cos A - (0.0341 + 0.0307M^) cos^ a] 


X (cos 


^ cos"^ 
60 


Vcos A + 1.992 


+ 0.0439 
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It is important to note the variation of grit drag increment with k Indicated by 
equations (A12) and (A13) when all other variables are held,constant. For k/6*j^s 5, 
equation (A12) indicates ACj^ varies as k^; whereas for k/6*j^ > 5, equation (A13) 
indicates ACj^ varies linearly with k. This discontinuity in the variation of AC^ 
with k has been observed from experimental data obtained by other investigators. 
Figure 12 of reference 2 provides a good example. The discontinuity in slope of the 
reference 2 data occurs at k » 56*^, and for k<56*j^, ACj^ varies as k^ which is 
consistent with predicted trends of the present method. 

Applying equations (A12) and (A13) to the model used in the present tests where 

A = 55° 

S = 0.2045 m2 

I = 243.7 cm 

results in the following simplified equations 

w2 

ACq = (0.00265 - O.OOOIM 

® k 
and 

ACj^ = (0.0133 - 0.0005Moo)k 
where the units of both k and 6*^^ are cm. 




> 5 


(A14) 


(A15) 
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(a) Model geometry; all dimensions in cm. 
Figure 1.- Model details. 
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(b) Photograph of model. 
Figure 1.- Concluded. 
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Figure 3.- Effect of unit Reynolds number on transition location and transition 
Reynolds number, k = 0 cm, a = 0°. 
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Figure 4. - Effect of angle of attack on transition location, k = 0 cm. 
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Figure 7.- Comparison of effective roughness Reynolds numbers 
from present tests with existing flat -plate data. T , = T„„ : 
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a = 0^; k ^ 6^. 
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Figure 10,- Effect of k on drag at zero lift, 
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Figure 10.- Concluded. 
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(a) Concluded. 
Figure 13.- Continued 
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Figure 14. - Comparison of calculated and measured drag coefficients for circular cylinders used as 

boundary -layer trips. 






Figure 15.- Comparison of calculated and measured drag 
coefficients for isolated cylinder in presence of 
laminar boundary layer, k = w. 
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